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ABSTRACT 

Recent analysis of five gamma-ray burst (GRB) afterglow spectra reveal the absence of molecular 
hydrogen absorption lines, a surprising result in light of their large neutral hydrogen column densities 
and the detection of H2 in similar, more local star-forming regions like 30 Doradus in the Large Mag- 
ellanic Cloud (LMC). Observational evidence further indicates that the bulk of the neutral hydrogen 
column in these sight lines lies more than 100 pc from the GRB progenitor and that H2 was likely 
absent prior to the burst, suggesting that direct flux from the star, FUV background fields, or both 
suppressed its formation. We present one-dimensional radiation hydrodynamical models of GRB host 
galaxy environments, including self-consistent radiative transfer of both ionizing and Lyman- Werner 
(LW) photons, nine-species primordial chemistry with dust formation of H2, and dust extinction of 
UV photons. We find that a single GRB progenitor is sufficient to ionize neutral hydrogen to distances 
of 50 - 100 pc but that a galactic LW background is required to dissociate molecular hydrogen in the 
ambient ISM. Intensities of 0.1 - 100 times the Galactic mean are necessary to destroy H2 in the cloud, 
depending on its density and metallicity. The minimum radii at which neutral hydrogen will be found 
in afterglow spectra is insensitive to the mass of the progenitor or the initial mass function (IMF) of 
its cluster, if present. 

Subject headings: cosmology: theory — galaxies: star clusters — gamma rays: bursts — H II regions — 
ISM: clouds — radiative transfer 



1. INTRODUCTION 

Long-duration GRBs are extremely energetic events 
produced during the death of some massive stars. Be- 
cause massive stars have such short lifetimes {t < 10^ 
yr), the GRB events are identified exclusive ly with star- 
forming galaxies (e.o. lLe Floc'h et al. II2003I ) and are ex- 
pected to occur within the star-forming regions of these 
galaxies. In addition to the burst of gamma-ray radia- 
tion, approximately half of high z GRBs emit large fluxes 
of rest-frame UV radiation observed as bright optical af- 
terglows at Earth. These afterglows are believed to re- 
sult from the deceleration of a relativistic jet upon its 
encounter with gas at r ~ 10^^ cm from the GRB pro- 
genitor. The optical afterglows have peak brightnesses 
which rival (and sometimes exceed) the brightest quasars 
at high z. 

In response to the discovery of optical afterglows, as- 
tronomers have initiated programs to obtain photome- 
try and spectroscopy at early times (t < I hi) before 
these transient events fade below detection. Imprinted 
on the power-law spectrum of a GRB afterglow are the 
absorption-line signatures of all the neutral and partially 
ionized gas between the GRB and Earth. In principle, 
this includes the circumstellar medium (GSM) associated 
with the GRB progenitor, gas comprising the molecular 
cloud and/or H II region related to the local star- forming 
region, the 'ambient' interstellar medium (ISM) which 
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surrounds the star-forming region, baryons in the galac- 
tic halo of the GRB host, and the so-called intergalactic 
medium (IGM) that permeates the universe. 

To date, astronomers have struggled to identify signa- 
tures of the expected CSM, most likely because the GRB 
afterglow it self ionizes this gas prior to spectroscopi c ob- 
servations ([ChinTIaD [2Q03; l^tioZiS^D [IQO^ A 
typical GRB afterglow that is bright enough for follow- 
up spectroscopy emits ~ 10^" ionizing photons during 
the first few minutes after onset (assuming isotropic 
emission). This implies ^ 10^° photons cm~^ at 10 pc, 
and previous calculations have shown the GRB ionizes 
and dissociates the majority o f H and H2 to this ra- 
diusJWaxman & Drami l200d : [P erna fc Lazzatil l2002t 
iDraine fc Hao 2002). At distances greater than 10 pc, 
however, the GRB afterglow should have minimal effect 
on the ionization state of the star-forming region and 
surrounding ISM. 

The gas observed within the host galaxy of GRBs is 
marked by large neutra l hydrogen column density, rifji 
(jJakobsson et all |2006( ) . Many of the sightlines show 
A^(H2) > lO^^cm"^, surface densities comparable to 
those measured for molecular clouds. A recent sur- 
vey for H2 along a subset of these sightlines, however, 
has demonstrated that the molecular fraction /(H2) = 
2A^(H2)/[2A^(H2) + A^(HI)] is ge nerally very low in the 
gas surrounding GRBs, < IQ-^ '^ (jTumiinson et al.ll2007| ) 
(hereafter TET07). The authors draw two conclusions 
based on this observational result: (i) the molecular gas 
associated with the GRB star-forming region was pho- 
todissociated prior to the GRB event; (ii) the formation 
of H2 in the ambient ISM is surpressed by an enhanced 
far-UV radiation field. However, we note a tentative de- 
tection of H2 has been made tow ard GRB 060206 with 

/(H2) 3.5 (|Fvnbo et al.l[2006l ). TET07 regarded this 

result as an upper limit because they consider the iden- 
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tifications of the two possible H2 absorption lines to be 
ambiguous. Even if we allow this tentative discovery, the 
absence of molecular hydrogen in other GRB afterglows, 
to very low detection limits, requires explanation. 

There are several lines of argument which suggest the 
gas observed in GRB sight lines arises primarily in the 
ambient ISM of the host galaxy. First, many afterglow 
spectra exhibit strong Mg I absorption at velocities iden- 
tical to the (dominant) low-ion transitions of S II, Zn II, 
Mg II, etc. Although the GRB afterglow will photoionize 
hydrogen and these low-ion species to only ^ 10 pc from 
the GRB. lProchaska. Chen fc BloomI (|2006[ ) argued that 
atomic Mg should b e ionized to beyond 50 pc (see also 
iMirabal e^ "ani2002[ ). The observation of strong Mg I 
absorption, therefore, places the majority of gas at dis- 
tances that probably lie beyond the local star-forming 
region. Second, the gas along GRB sight lines also tends 
to exhibit tra nsitions from fine-structure level s of Si"*", 
0°, and Fc+. (iProchaska. Chen fc Bloomll2006l ) demon- 
strated that this gas is excited by UV photons of the 
GRB afterglow, which sets an upper limit on its distance 
from the event of 1 kpc. Third, UV pumping lends 
to line variability as the afterglow fades which allows a 
more precise estimation of the distance of the gas from 
the afterglow. Analysis of GRB 020813 and GRB 060418 
indica tes distances of ^ 100 p c and 1.7 kpc respec- 
tively (jPessauges-Zavadskv et al . 2006; Vrccswiik et al. 
120071 ). Finally, as pointed out by .Watson et al. . (,20071 
independent x-ray and optical evidence also suggests that 
the neutral columns are at some distance from the GRB 
but local to the source galaxy. To summarize, the obser- 
vation and analysis of Mg I and fine-structure absorption 
indicates the majority of gas observed in the afterglow 
spectrum is located at ^ 100 pc to 1 kpc. 

Taken together, these results inspire the following ques- 
tions: first, why is there so little molecular or even neu- 
tral gas within ~ 100 pc of the GRB? Does this imply an 
extensive H II and photodissociation region (PDR) , and 
if so what is the source of this radiation field? If a sub- 
stantial far-UV flux is suppressing II2 formation in the 
ambient ISM prior to the burst, is the star- formation re- 
gion associated with the GRB responsible or were other 
sources required? 

We present numerical simulations of the photoioniza- 
tion and photodissociation of a molecular cloud by a mas- 
sive star cluster hosting a GRB progenitor. We survey 
the typical radii of the H II region and central PDR for 
reasonable collections of stars and evaluate the role of 
the cluster and galactic UV background in suppressing 
H2 formation in the outer regions of the cloud. Our pa- 
per is organized as follows: in § 2 we describe the numer- 
ical algorithm applied to the GRB progenitor H II region 
and PDR. The evolution of the central H II region and 
II2 abundances in the outer cloud are examined for a 
fiducial case in § 3. The photoionization of the molecular 
cloud core is computed for a variety of central densities 
and fiuxes in § 4, and the minimum LW backgrounds re- 
quired to destroy H2 in the cloud are tabulated in §5. 
We summarize in § 6. 

2. NUMERICAL ALGORITHM 

We evolve the molecular cloud hosting the GRB with 
ZEUS-MP (iWhalen fc NormanI [20061 l2007ah . a mas- 
sively parallel Eulerian reactive flow hydrocode with self- 



consistent m ultifrequency photon cons erving UV radia- 
tive transfer (jWhalen fc Normanl2007bl ) and nine-species 
primordial chemistry. ZEUS-MP can transport photons 
from a point source centered in a spherical grid or plane 
waves along the x- or z-axes of cartesian or cylindrical 
boxes. Updates to the gas energy equation include cool- 
ing due to electron coUisional excitation and ionization 
of H, recombinations, bremmstrahlung, inverse Comp- 
ton scattering and ro-vibrationa l emission by molecu- 
lar hydrogen (Galli & Palla 199 q). We also in corporate 
metal line cooling (D algarno fc McCravlllQTl) scaled to 
the metallicity of the gas. 

Our reaction network is the 9 -speci es pri mordial gas 
phase model of lAnninos et al. I (|1997D and lAbel et oTI 
(|1997f ). We augmented the gas phase chemistry with 
dust c atalysis of II2 according to the prescription of lJural 
(fT97l 

^ ZRnnnu, (1) 

where R is estimated empirically to be ~ 1 - 3 x 10~^^ 
cm~^ in the solar neigborhood and Z is the metallicity 
in units of Zq. In doing so, we make the tacit assump- 
tion that grain formation of II2 is independent of tem- 
perature. Although this approximation is widely applied 
in PDR models (and is reasonable for Tdust ^ 50 K), 
it should be recognized that R probably falls with ris- 
ing dust temperatures beca use the 'sticking factor' of H 
atoms on grains decrea ses (jHollenbach fc McKed 119791 
ICazaux fc Tielenj|2002[ ). The temperature dependence 
of R remains poorly constrained because the grain struc- 
ture is unknown. This is not an issue within H II regions 
since dust there will be photoevaporated; however, be- 
cause the rest of the cloud could be heated above 100 K 
by a galactic background, our results for II2 production 
there should be considered an upper limit. 

2.1. Photon-Conserving Dust Extinction 

We modified our multifrequency UV transport to ac- 
comodate dust attenuation using the extinction curves of 
[Pei (1992) in a manner that preserves photon conserva- 
tion. In photon-conserving radiative transfer, the total 
number of absorptions in a grid zone per second is the 
number of photons entering the zone per second minus 
those exiting: 

iVabs cx 1 - e-^ (2) 

where 

I 

T = ^ (JtUiAr (3) 

i=l 

is the total optical depth of the cell (the sum is over all 
removal processes). In eq[31 Ui is the number density of 
the species absorbing the photon in the i*'' reaction, ai 
is the interaction cross section of the reaction, and Ar 
is the zone length in the direction of propagation. We 
extract the number of absorptions Ni due to a particular 
reaction from TVabs by 

N. = N^.s^0^y (4) 
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where the ti are the optical depths to aU the processes 
and Ti is that associated with the reaction. 

We incorporate dust absorptions into this procedure 
by adding td, the local optical depth due to dust in the 
cell, to the sum in eq[3l evaluating A'^abs with eqO and 
adding 1 - td to the sum in the denominator of eqlH 



1 - e- 



abs 



(l-e-B) + ^^(l_e--0' 



(5) 



The effect of dust absorptions is to reduce the number 
of all other absorptions in the mesh zone. The radiative 
rate coefficients used in the reaction network are then 
easily derived from the Nf. 

h = (6) 

2.2. Extinction Curves 



We obtain td for a mesh zone from m (UMl) as fol- 
lows: 

rD(A) = rBe(A), (7) 

where tb is the extinction optical depth in the B band 
and Td includes both absorption and scattering. Dust 
and gas are not selectively transported in our model so 
Tb can be defined in terms of an observationally derived 
dust-to-gas ratio k 



where 



(8) 



(9) 



is the total column density of hydrogen. The scale factor 
^ is given by 



- (A/A,)"' + (Ai/A)". +5/ 



(10) 



We adopt the dust-to-gas ratio from Table 2 and pa, - 
rameters a^, 6i, n^, and Ai from Table 4 of iPeil ()1992f ) 
for the Small Magellanic Cloud (SMC) for our opacity 
model, appropriate for the meta llicities that characteriz e 
the host galaxies of GRBs {e.g. iProchaska et al. ]l2007l ). 
We initialize free electron fractions in the simulations by 
assuming each of the metal species in the Dalgarno & 
McCray curves is singly-ionized by a background UV ra- 
diation field. 

2.3. Multifrequency UV Transfer/ Uniform LW 
Background 

We partition photon emission rates by energy bins ac- 
cording to the blackbody spectra of the stars in the clus- 
ter. Good convergence in frequency is found with 40 
uniformly-spaced bins in energy from 0.755 eV to 13.6 
eV and 80 logarithmically spaced bins between 13.6 eV 
and 90 eV. Flux in the lower energy range cannot ionize 
H or He but does drive a host of chemical reactions that 
regu late H~ and H2 for r nation and are listed in Table 
1 of IWhalen fc NormanI (|2007br ). H2 photodissociation 
rates are computed with the self-shielding functions of 



iDraine fc Bertoldil (Il996l) modified for thermal doppler 
broadening (jWhalen fc NormanI [2007c? ) to address gas 
flows in an approximate manner. Multiple stars can be 
stationed at the origin but they are not spatially resolved 
and must all be of the same mass (and blackbody spec- 
trum). Either constant luminosities or time-dependent 
fluxes derived from stellar evolution models can be used. 
Stellar spectra can be updated on Courant times without 
serious loss of efficiency even though they usually evolve 
on much longer timescales. 

We approximate galactic FUV dissociating back- 
grounds with uniform fluxes that permeate the cloud, 
which are specified in multiples of J21 = 1 x 10~^^ erg 
cm"2 fj2-i str-i s"i (jMachacek et al.ll200lD . They are 
implemented by a simple addition to the photodissocia- 
tion rate coefficient in the reaction network 



1.13 X 10^ 47rJ2i. 



(11) 



The Milky Way FUV background intensity (2 x 10"^ 
photons cm~^ s~^ Hz"-'^, centered at hv = 12.87 eV) is 
related to J21 by /mw — 32.82 J21. 

3. H II REGION/PDR DYNAMICS 

We examine first the breakout of UV radiation from 
a cluster of massive stars in which a GRB progenitor 
resides. The stars are situated in a molecular cloud core 
of density surrounded by a uniform envelope density 
rici. The density profile of the core is represented by a 
radial power law: 



n{r) 



ric r <rc 

n-cir/r^y^ r > r^, 



In all the models in this paper we set r^ = 
0.2 pc and ut — 2, c onsistent with observ a tional 
surveys of cloud cor es JArauilla fc GoldsmithI [l985l : 



iGregorio Hetem et al.lll988D . We smoothly ioin this pro- 
file to the uniform density Ud assumed for the rest of the 
cloud at the radius at which they become equal. This 
second radius varies with Uc and Ud and constitutes the 
outer boundary of the core, since it is the distance at 
which the core becomes indistinguishable from the cloud. 
The size of the cloud itself is just the outer boundary of 
the grid, as discussed below. Here, n, Uc, and are 
total particle densities, partitioned by the mass fractions 
of the species assumed initially to be on the grid, 0.76 
H, 0.24 He in the models of this section. We set this 
gas in hydrostatic equilibrium by computing the gravi- 
tational potential whose forces exactly balance pressure 
forces everywhere on the grid. Gravitational force up- 
dates computed from this potential are applied to the 
gas velocities at every hydrodynamical time step in a 
run but the potential itself is held fixed. The self-gravity 
of the gas is not included in our models because it is neg- 
ligible in comparison to the pressure forces in the ionized 
gas. 

A source of ionizing flux was placed at the center of a 
one-dimensional spherical mesh, whose inner and outer 
boundaries were 0.1 pc and 1000 pc. Reflecting and 
outflow conditions were assigned to the inner and outer 
boundaries, respectively. The grid was discretized into 
200 logarithmically-spaced bins deflned by 
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r(pc) r(pc) r(pc) r (pc) 

Fig. 1. — Panels (a) - (d): density, temperature, ionization fraction, and velocity profiles of the photoionized GRB host cloud. From left 
to right in each panel, the curves are at t = 38.7 kyr, 773 kyr, 1.55 Myr, 3.1 Myr, and 3.87 Myr, respectively. Panels (e) - (h): final ionized 
fraction, density, temperature, and H2 fraction profiles at 3.87 Myr for the five physics test models. Black: model 1; blue: model 2; red: 
model 3; green: model 4; magenta: model 5. In panels (e) - (g), models 1-3 largely overlap at the larger radii while models 4 and 5 (in 
which dust extinction is included) overlap at smaller radii. In panel (h), H2 fractions are below the scale of the plot. 



— P ^1 — ^ outer ^inncri (-^^) 

n 

where n is the number of radial zones. We apphed a grid 
ratio (3 = 1.043, which resolves the flat central core with 
10 zones. 

3.1. The H II Region 

Ten 40 Mq stars were centered in a molecular cloud 
with lie = 100 cm~'^ and rid = 5 cm~^, with a neutral 
H column of ~ 1.5 x 10^^ cm^, consistent with those of 
the GRB sight lines. The ionizing flux of the cluster will 
reduce this neutral column but never by more than 15 % 
because the H II region rarely grows beyond 150 pc, as 
we will show below. The temperature and initial H2 frac- 
tion of the gas is 100 K and 2 x 10"'', respectively, and 
the UV source is maintained for the 3.87 Myr lifetime of 
a 40 Mq star. For simplicity, we assume a uniform ion- 
izing luminosity, effective surface temperature, and total 
luminosity of 2.469 x 10"*^ photons s"\ 7.94 x 10"* K, 
and 2.63 x 10^ L©, respectively, for each star according 
to the zero-metallicity models of Schacrcr (2002) . We in- 
clude a LW background equal to that of the Milky Way, 
a metallicity of 0.1 Zq, and both production of II2 and 
extinction of photons by dust. 

The ionizing UV photons of the cluster initially propa- 
gate into the cloud behind an ionization front (or I- front), 
whose thickness is only ~ 20 photon mean free paths 
(mfp) in the neutral gas. The neutral fraction abruptly 
transitions from 10~^ to 1 across this layer, which at first 
expands supersonically with respect to the ambient gas. 
This type of I- front is known as an R-type front, and 
its passage leaves the gas ionized and hot, but other- 
wise undisturbed. However, the I-front eventually slows 
down due to geometric dilution of the photons reaching 



the front and recombinations in the ionized gas, or H 
II region. Recombinations to the ground state absorb 
a photon from the cluster but release another ionizing 
photon, with no net effect on the advance of the front. 
Recombinations to excited states absorb photons from 
the cluster without returning any ionizing photons to the 
I-front and slow its expansion. In a static medium, when 
the number of recombinations to excited states in the 
entire H II region equals the number exiting the cluster, 
the I-front has reached the Stromgren radius and can 
advance no further. 

In a real gas, a pressure wave builds in the H II re- 
gion. As the pressure wave accelerates and the front de- 
celerates, the wave overtakes the front at approximately 
the Stromgren radius, steepens into a shock, and breaks 
through the front. The pressure of the ionized gas drives 
the neutral shock out into the cloud. As the shock grows 
in radius, it sweeps up gas in the cloud, which accumu- 
lates in a dense shell that detaches the shock from the 
I-front. This kind of I-front is known as a D-type front, 
and it advances supersonically with respect to the cloud, 
but subsonically with respect to its own ionized gas. 

A good approximation for the initial Stromgren radius 
of the H II region (which accounts only for hydrogen 
ionizations) is 

where riph is the emission rate of ionizing photons in s^^ , 
q;b is the case B recombination coefficient of H, and nn 
is the number density of H. The final radius of the H 
II region (after coming to pressure equilibrium with its 



5 



surroundings) is 



where and Tq are the ionized gas and ambient tem- 
peratures, respectively. Given our model parameters, rg 
~ 12 - 90 pc, depending on whether core or cloud densi- 
ties are used, and r/ is 85 times greater. Thus, the H 
II region never comes into pressure equilibrium in the life 
of the association and is still expanding when the GRB 
explodes. As we later show, dust extinction reduces r^ 
by a third in our runs. We adopt a coohng cutoff of 
1000 K, defined to be the gas temperature below which 
cooling curves are not applied to the gas energy updates. 
In reality, the Dalgarno & McCray curves would reduce 
the temperature of the shocked shell driven by the ion- 
ization front (I-front) to below 100 K, collapsing it to 
higher densities and narrower widths than in our mod- 
els. However, this results in less forgiving Courant times 
and alters neither the escape of UV nor the dynamics of 
the H II region, so we proceed with the higher cutoff. 

Density, temperature, ionization fraction, and velocity 
profiles are shown in Figs[T^-d at t = 38.7 kyr, 773 kyr, 
1.55 Myr, 3.1 Myr, and 3.87 Myr. The I-front at first 
propagates so rapidly that no hydrodynamical response 
of the gas can overtake it. This is an R-type front, and 
its supersonic nature is evident at 38.7 kyr: the gas is 
barely disturbed and the growing velocity wave lags the 
front. As the I-front expands it slows down due to both 
geometrical dilution of radiation and recombinations in 
the ionized gas that reduce the number of ionizing pho- 
tons reaching the front. The pressure wave building up 
in the hot interior catches up to the front, pushes past 
it, and steepens into a strong shock. The I-front is now 
D-type, and as it continues to grow it sweeps up shocked 
neutral gas in a dense shell that detaches the shock from 
the front. In this model, the I-front executes a transition 
from R-type to D-type at 460 kyr. The relatively large 
postfront gas temperatures of ~ 35,000 K are due to the 
ionization of He by the hard UV spectrum of the star 
cluster. In the large densities of the core, the spread in 
mean free paths of the UV photons is fairly narrow so the 
I-front is barely resolved. As it reaches the lower intra- 
cloud densities the front broadens to approximately 5 pc. 
After becoming D-type, the I-front drives the shock into 
the neutral medium at 13 km s~^. As the H II region 
expands, its densities fall from 100 cm~'^ to less than 1 
cm~^, so the GRB explodes in a comparatively diffuse 
environment. By 3.1 Myr the shock separates from the 
I-front, as shown in the temperature profile. The densi- 
ties in the neutral shell snowplowed by the front rise to 
10 cm~'^. 

At the cooling cutoff temperature the shocked shell 
remains fairly thick, approximately 20 pc. Even so, in 
three dimensions the I-front would be prone to violent 
dynamical instabilities me diated by radiative cool ing in 
the shell by metal lines ((Whalen fc Normanll2007af l. The 
cooling fractures the shell, allowing radiation to jet for- 
ward through the fissures in prominent fingers of ion- 
ized gas much larger than the original perturbations. 
These i nstabilities would form eve n if no metals were 
present (jWhalen &: NormanI l2007"bl ) because the I-front 
forms H2 at the base of the shocked neutral shell. As 



TABLE 1 
Five GRB Host Cloud Models 



model 


Xe 


dust H2 


dust extinction 


LW Background 


1 


no 


no 


no 


no 


2 


yes 


no 


no 


no 


3 


yes 


yes 


no 


no 


4 


yes 


yes 


yes 


no 


5 


yes 


yes 


yes 


yes 



we show in the model 5 curve in Fig [TJi, molecular hy- 
drogen formed via the H^ channel in the partially ion- 
ized outer layers of the broadened I-front is not entirely 
destroye d by Lyman- Werner flux fr om the cluster (sec- 
tion 2 of IWhalen fc NormanI l2007bl ) . The 2 x lO'^ H2 
abundances, although not as efficient as metals, still suf- 
ficiently cool the sh ell for instabilities to erup t in the 
Tfront (see Fig 5 of IWhalen fc Normanl[2007H l. These 
features are relevant to the post-afterglow evolution of 
the GRB because they form clumps that would be ad- 
vected through the shock front. However, it is unlikely 
that these clumps would harbor enough H2 for detection 
along sight lines. 

3.2. The Photodissociation Region 

We performed five runs in which elements of our 
physics model were sequentially adopted in order to eval- 
uate their individual impact on H II region dynamics, 
the central PDR, and equilibrium H2 abundances in the 
outer cloud. In the first run no initial electron fraction 
(Xe), dust production of H2, dust extinction of photons, 
or Lyman- Werner background radiation were present. 
The H II region and PDR would be largest in this model 
because no dust absorbs photons and there are no free 
electrons to restore H2 lost to dissociation by the cluster. 
In the second model, an initial Xe consistent with singly- 
ionized metal species was included, so H2 levels beyond 
the H II region attain an equilibrium set by photodisso- 
ciation by central LW photons and formation by the H^ 
and H^ channels in the gas phase: 

H e- ^ H- + 7 H- + H ^ H2 + e" 
H++H ^ H++7 H++H ^ H2+H+ 

The H~ channel is dominant in these density and tem- 
perature regimes and will restrict the final size of the 
PDR by replenishing H2 lost to LW photons from the 
center. The activation of H2 formation by dust grains 
in model 3 further reduces the PDR. Neither of these 
processes affect the size of the H II region, but activa- 
tion of UV absorption by dust in model 4 will slow the 
front. In model 5 we add a uniform LW fiux equal to the 
Milky Way background that dissociates H2 everywhere 
in the cloud. It is the fiducial case discussed above. We 
summarize the five models in Table [TJ 

Ionized fractions, densities, temperatures and H2 
abundances for the cloud at 3.87 Myr appear in Figs[Tt-h. 
The fiow profiles clearly fall into two classes: those that 
incorporate dust extinction and those that do not. We 
find that dust attenuation at abundances similar to the 
SMC reduce the radius of the H II region by 30 % at the 
end of the main sequence lifetime (MSL) of the 40 Mq 
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Fig. 2. — Panel (a); emission rate of all ionizing photons as a function of time. Panel (b); radii of the cluster H II regions at the time of 
the gamma-ray burst. Filled symbols: 25 Mq cluster; unfilled symbols: 40 Mq cluster. Triangles: nd = 5 cm~^; squares: rid = 10 cm~^; 
diamonds: = 50 cm~^. 

star, from 120 pc to 85 pc. Lower photoionization rates 
drop postfront gas temperatures by a few thousand K, 
which in turn shghtly reduces the shock velocity. Aside 
from this there are few differences in the structure of the 
flow curves of the two groups: shocked shells and I- fronts 
have nearly the same width at 3.87 Myr. 

On the other hand, each additional pliysical process 
significantly alters H2 fractions in the cloud. The dashed 
curve represents model 1, in which no free electrons or 
dust form H2, and it peaks at 125 pc (being overlaid 
by the model 2 and 3 curves) but then drops to very 
low levels beyond. The peak is the molecular hydrogen 
manufactured in the front; it cannot shield the initial 2 
X 10~^ H2 fractions in the outer region of the cloud from 
LW flux from the cluster, so they are mostly destroyed. 
Molecular hydrogen within the H II region is coUisionally 
dissociated by free electrons to extremely low levels. The 
small free electron fraction (~ 3.5 x 10~^) maintains a 
little of the original H2 beyond the H II region at 500 
pc in model 2 as shown by the dotted curve, but the 
entire cloud is still essentially dissociated by the cluster. 
Activation of II2 production by dust shields the outer 
cloud from LW flux, and molecular hydrogen abundances 
there rise by an order of magnitude above their initial 
values beyond 500 pc. The PDR in this model has a 
radius of 250 pc. 

In model 4, dust attenuation of photons absorbs all 
central LW flux within 5 pc of the H II region, so H2 
fractions beyond rise by a factor of ten in 3.87 Myr, to 2 
X 10""^. Molecular hydrogen abundances in the 1-front 
are 100 times greater than in models 1, 2, and 3 due to 
FUV absorption by dust in the H II region. In reality, 
this peak in H2 would subside over time as the dust is 
evaporated by ionized gas, so the PDR would expand 
slightly more than 5 pc beyond the H II region. The 
uniform LW background activated in model 5 reduces 



II2 fractions in the 1-front peak to 2 x 10^^ and to 2 
X 10^^ in the rest of the cloud, lowering H2 mass frac- 
tions below 10~^'^, the detection limit in the GRB sight 
lines. All four processes are key to constructing accurate 
molecular hydrogen abundances in the GRB host cloud. 
Models 4 and 5 clearly demonstrate that H2 abundances 
in the cloud are set by the LW background, since for re- 
alistic models of the cluster the central PDR is confined 
to within a few pc of the H II region by dust absorption 
(rather than by H2 formation in cither the gas phase or 
on dust grains). Our models demonstrate that the sizes 
of the H II regions, and therefore both the column den- 
sity of neutral H and the minimum radii at which it is 
observed, are decoupled from what sets H2 abundances 
in the cloud. In the next section we examine the mini- 
mum distances from the GRB at which neutral H would 
be detected for a variety of clusters and cloud core densi- 
ties. In the section thereafter, we evaluate the minimum 
galactic UV backgrounds required to drive H2 below de- 
tection limits. 

4. THE GRB HOST CLOUD H II REGION 
4.1. Clusters 

For what range of central densities and stellar masses 
is the molecular cloud core completely ionized, and what 
final H II region radii result? We consider two clusters 
of ten stars, one composed of 25 Mq stars and the other 
with 40 Mq stars. For rid = 5, 10, and 50 cm~^ we 
consider central core densities of 10, lO^, 10'\ 10"*, and 
10"'' cm^'', a total of 30 models. The ric are consistent 
with molecular cores observed in the Galaxy, and rid = 
5-10 cm~'^ are typical of clouds in the Milky Way and 
yield the neutral hydrogen columns of the GRB spectra 
given our mesh boundaries. We include rid = 50 cm~^ as 
a reasonable upper limit in case average cloud densities 
rise with redshift, 
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For consistency with the metaUicity of the GRB host 
galaxy, we adopt the time-dependent spectra of 0.1 Zq 
25 and 40 M© stars in this parameter survey. At a 
given mass, a primordial star is more compact and hot- 
ter than a 0.1 Zq star for two reasons. First, Population 
III stars derive their nuclear energy primarily from in- 
efficient proton-proton burning, with only a small con- 
tribution from the CNO cycle due to small amounts 
of C that are built up in the pre-rnain se quence phase 
(|Castellani et all ITosl lEl Eid et allflOSl . Since at a 
given temperature this process has lower energy genera- 
tion rates, the core must rise to higher temperatures to 
support the star from gravitational collapse. This, to- 
gether with the fact that there are no metals to absorb 
photon momentum in the photosphere and expand and 
cool it, makes primordial stars consi derably hotter than 
meta l-enriched stars of equal mass (jTumlinson fc Shulll 
I2OOOD . The time-dependent spectral fluxes for these stars 
were computed using the KEPLER code (jWeaver et al.l 
assuming a time-dependent black body spectrum 
at the effective temperature of the star. The result is 
shown in Fig [2^. At end of central hydrogen burning a 
very brief rise in the ionizing flux due to contraction of 
the star is observed, but then the stars become giants and 
the ionizing flux drops quickly to levels that are essen- 
tially negligible when compared to those when the star 
is on the main sequence. We implement these spectra 
in the code as photon emission rates that are binned by 
energy, as described earlier. Notice that the to tal emis- 
sion r ates of both stars are lower than those of iSchaereil 
1200^ because they are somewhat cooler, and that their 
luminosities can vary by 30% over their lifetimes. The 
0.1 Zq stars also have longer lives: 5.65 Myr vs 3.86 Myr 
for 40 Mq and 8.32 Myr vs 6.46 Myr for 25 Mq. The co- 
ordinate grid and all other model parameters were those 
used in section [3TT] 

H II region radii at the time of the GRB for all 30 runs 
are shown in FigOs. Several trends are apparent. First, 
the cloud density rid determines the H II region radius 
when ric is small because the front quickly ionizes the 
core. When central densities are greater than 10"* cm'^, 
the front requires a significant fraction of the lifetime of 
the star to break out of the core so Uc begins to also 
regulate the final radius. At any core density, the H II 
region shrinks with larger rid due to its greater inertia, 
as expected. Interestingly, for a given ric and Ud the less 
massive cluster has the larger H II region. Overall, the 
H II region radius at the time of the GRB varies little: 
from 41.1 to 100.7 pc over four orders of magnitude in 
He, an order of magnitude in rid, and a factor of five in 
luminosity. 

This somewhat surprising lack of variation in radius is 
due to two factors. First, if the source spectrum can ion- 
ize helium then postfront gas temperatures are ^ 25,000 
K, varying little with the number of stars in the clus- 
ter and only somewhat with their surface temperatures. 
This confines the sound speed Cg of the ionized gas, and 
therefore the expansion rate of the D-type front, to a 
fairly narrow range because of the weak dependence of 
Cs on T. Since the I-front will become D-type in a small 
fraction of the life of the star in realistic clusters because 
of dust, the velocity of the front is determined by the 
sound speed of the gas rather than the size or luminosity 
of the cluster. Dust augments the role of recombinations 
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Fig. 3. — Radii of H II regions with constant luminosities (unfilled 
symbols) and time- varying spectra (filled symbols). Triangles: 40 



Mq stars with n^i 
_3 



10 cm ■^i squares; 25 Mq stars with n^i = 5 



in the ionized gas by absorbing photons before they can 
reach the I-front and advance it, hastening its transition 
to D-type. 

Second, the time for which the H II region has grown 
when it is illuminated by the afterglow is, to first order, 
the lifetime of the progenitor because stars in the clus- 
ter form at approximately the same time. Since main 
sequence lifetimes for likely progenitors vary by no more 
than a factor of two, we expect neutral H in general to 
be observed 50 - 100 pc from the GRB. We note that 
since the final radius is determined by both the life of 
the cluster and the sound speed of the ionized gas, the 
25 Mq cluster has the larger H II region simply because 
its stars live longer. 

A simple estimate of the size of the H II region at the 
time of the burst follows from 



(15) 



For completely ionized gas (/i — 1.70) at 25000 K, Cg ~ 
14.3 km s-^ With MSL ranging from 3.87 to 6.46 Myr, 
the final radii of the H II regions would fall between 56 
and 94.3 pc, similar to those found by ZEUS-MP. 

Comparison of the sizes of the H II regions in models 3 
and 4 of section 13 . II suggests that the H II regions of clus- 
ters in higher redshift galaxies as a rule may have been 
larger than those in the Galaxy today because there was 
less dust to attenuate photons. We neglect the evapora- 
tion of dust within the H II region and therefore some- 
what underestimate the flux that reaches the ionization 
front, so these are lower limits for their radii. However, 
models 3 and 4 in section [XT] show that this effect is at 
most 30 % (and probably less, because dust on both sides 
of the I-front is not immediately destroyed and will still 
retard it). 
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Fig. 4. — Radii of the progenitor H II regions at the time of the 
gamma-ray burst. Filled symbols: 25 Mq star; unfilled symbols: 
40 Mq star. Triangles: n^i = 5 cm~^; squares: rid = 10 cm~'^; 
diamonds: n^i = 50 cm~'^. Note that the sizes of the H II regions 
of both progenitors match in n^i = 5 cm~'^ for rtc < 10** cm~^ 

The Milky Way LW background assumed in all 30 runs 
efficiently dissociates any H2 in the higher densities be- 
tween the H II region and the outer boundary of the core. 
Beyond this boundary, molecular abundances assume the 
same low equilibrium values as in model 5 of section [3T2l 
regardless of 11^. Again, we emphasize that these are not 
the ultimate sizes of the H II regions because, far from 
being in pressure equilibrium with their environments, 
they continue to expand as long as there is sufficient flux 
to maintain ionization. These are simply their radii at 
the time of illumination by the GRB. 

4.2. Stellar Spectra 

We have shown that the masses of the stars in the 
cluster have little effect on the final sizes of GRB host 
H II regions. Now, for completeness, we examine the 
H II regions of both zero-metallicity and 0.1 Zq stel- 
lar stars. This range in metallicity effectively brackets 
that of the host galaxies of the observed afterglows: 5.5 
X 10"'^ - 0.234. The variation in I-front evolution over 
these metallicities will probably be less than that over 
the stellar masses chosen for this study, since tempera- 
tures and luminosities change less with composition than 
with mass. In doing so, we also address in an approx- 
imate way the consequences of assuming a black body 
radiation spectrum without correcting for the stellar at- 
mosphere. Changes in the flux due to the presence of 
an atmosphere are less severe than the difference in flux 
between 25 and 40 M© stars. 

In Fig [3] we compare the final H II region sizes for both 
metal-free and 0.1 Zq 25 and 40 Mq clusters. We take 
the emission rate of ionizing photons, surface tempera- 
ture, total luminosity, and main sequence lifetime of the 
25 Mq stars to be 7.583 x s'\ 7.0 8 x 10^ K, 7.7 6 
lO'' Lq, and 6.459 Myr, respectively (|Schaereij|20"ol ). 



The 25 Mq clusters are in cloud densities ric 



5 cm 



and the 40 Mq clusters are in rici = 10 cm"'^. It is clear 
that time variation in the luminosity alters the final sizes 
of the H II regions by at most 15 %. The 0.1 Zq clusters 
have slightly larger H II regions because of their longer 
lives. Since the sizes of the H II regions are determined 
by the lifetimes of the stars and the sound speed of the 
ionized gas, H II regions of clusters with metals should 
be larger than those of clusters without metals by a fac- 
tor equal to the ratio of their lifetimes. This is greater 
than what we found because the spectra of the 0.1 Zq 
clusters soften as their stars cool with age, lowering Cg 
and slowing the expansion of the I-front. 

4.3. Solitary Progenitors 

So far, we have shown that as long as there is enough 
UV flux to break out of a molecular cloud core and es- 
tablish an H II region, its final size is largely indepen- 
dent of total flux or spectral profile. Could then a single 
GRB progenitor create the H II regions implied by the 
minimum distances at which neutral H is found in the 
afterglow spectra? We repeat the models of section |4] 
but use just one star instead of 10. The sizes of the H II 
regions are shown in Fig|3) 

Interestingly, the H II regions are similar to those of 
the clusters, attesting to the fact that postfront gas tem- 
perature rather than fiux governs their sizes. They also 
follow the same trends: their radii fall off somewhat with 
increasing at a given cloud density and even more 
sharply with rid- Final radii are greater for smaller stars 
due to their longer lives. The H II region of the 25 Mq 
star is slightly more sensitive to the core density, partic- 
ularly when cloud densities are small. We conclude that 
solitary progenitors are as effective as clusters in forming 
H II regions in the GRB host galaxy. 

Our H II region radii are insensitive to Tc and realistic 
choices of density gradient uj. In section 3.1 we showed 
that even if H II regions expand in core densities (10^ 
cm~^) rather than in cloud densities (10 - 50 cm^'^), 
they reach Stromgren radii > 10 pc and final equilib- 
rium radii that are many times larger. It is clear that 
extending rc to 1 pc, which is unrealistically large for 
the centers of galactic molecular cloud cores, would have 
very little effect on final H II region radii. Furthermore, 
the I-front evolutio n is very similar for density gradients 
of 1.5 > cj > 2.5 (jFranco et al.lll990D . so there wiU be 
little variation in final H II region radii over reasonable 
ranges of core density profiles. 

4.4. Dust Extinction 

The simulations described so far all assume SMC dust 
extinction for UV transport. How much smaller would 
the H II regions be if more dust is present in the GRB 
host galaxy? We calculated the H II region of a singl e 25 
AIq star with the strongest attenuation curve from iPeil 
(1992), which is that of the Milky Way (MW). The SMC 
and MW curves bracket the dust extinction that could 
realistically be expected in the afterglow spectra. The 
star was placed in the largest core and cloud densities 
in our survey (10 and 10^ cm~^, respectively). We find 
that the final radius of the H II region does not change 
with the new attenuation curve. Since this choice of pa- 
rameters is the most hostile to the formation of an H II 
region, we conclude that none of our other models would 
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Fig. 5. — Minimum LW backgrounds required to suppress H2 
mass fractions below lO"'''^ in our models, in units of the Milky 
Way average background. Filled symbols: ZEUS-MP; unfilled sym- 
bols: analytical. Triangles: n^i = 1 cm""^; squares: n^i = 5 cm""^; 
diamonds: = 10 cm~'^; circles: the uniform LW backgrounds 
implied by eg llTl for the afterglows in TET07, assuming diameters 
of 10, 100, and 1000 pc for the GRB host cloud. The triple for each 
afterglow is connected by a vertical dotted line and is designated 
by its GRB. 

be affected over the range of extinction lifcely to be found 
in the galaxies. 

4.5. Winds 

In our models we have excluded the 500 - 1000 km 
line-driven winds of the stars. Winds would at least par- 
tially blow away the molecular envelope of the star, ex- 
cavating a cavity of diffuse coronal gas at several million 
K and plowing a strongly shocked (and completely ion- 
ized) gas shell out into the cloud. However, they would 
not alter the size of the H II regions at the time of the 
blast. While it is true that winds would disperse dense 
cores and facilitate breakout of UV into the cloud, there 
is little variation in final H II region diameters with Uc 
anyway. Once the H II region encloses the bubble, it will 
continue to expand in spite of the dense shell driven by 
the wind because it is so hot it cannot recombine and 
attenuate the UV flux reaching the I- front. For plausible 
cloud densities, the shell will not overtake the H II region 
(100 pc) prior to the burst. Furthermore, dust would be 
evaporated inside the bubble but survive beyond it, with 
little difference to molecular fractions in the rest of the 
cloud. Because the properties of the bubble are not ob- 
servationally constrained by current afterglow data and 
their inclusion would change neither the minimum radii 
at which neutral H is observed nor the H2 column den- 
sities in the sight lines, we do not include winds in this 
study. 

Nevertheless, two important points about wind blown 
bubbles merit further investigation. First, interactions 
between strong ionized flows in the H II region and the 
shell would lead to complex morphologies in the circum- 
stellar environment of the GRB that could be crucial to 



its evolution (e.g. iFrever et al.ll2003l ). Future observa- 
tions might unveil these geometries. Second, UV and x- 
ray photons from the hot wind- blown shock could modify 
the structure of the H II region ()Cho fc Ka ng 2008). The 
hard radiation broadens the I-front and crea tes a tem- 
perat ure front (T-front) precursor ahead of it (|Qiu et al.l 
I2OO7D . The front widens because of the spread in pho- 
ton mean free paths in the high-energy tail of the x-ray 
spectrum. The T-front is the layer of hot, nearly neutral 
gas formed when the most energetic photons in the spec- 
trum penetrate far into the gas upstream of the front, 
ionize it and create secondary electrons capable of other 
collisional ionizations. This minute but energetic flux 
can only sustain ionized fractions of 10~^ - 10~^ ahead 
of the front, but this is sufficient to raise their temper- 
atures to 10"* K, considerably hotter than the shock it 
precedes. X-rays can either discourage or promote H2 
formation in the outer layer of the H II region. The par- 
tial ionization and warm temperatur es there are ideal fo r 
catalysis of II2 via the H" channel ( iRicotti et al.l[200l[ l. 
but the high temperatures of the T-front may curtail its 
formation. 

5. PHOTODISSOCIATION OF THE CLOUD BY THE 
GALACTIC BACKGROUND 

We now explore the LW intensities necessary for II2 
destruction for a variety of column densities and metal- 
licities that are consistent with observations of the GRB 
progenitor cloud. As noted in section 12. 3i we do not 
perform direct transport of galactic UV photons in the 
outer regions of the cloud so we make no attempt to 
model the rise of either LW backgrounds or II2 fractions 
there. Molecular hydrogen, dust, and dissociating pho- 
tons all come to equilibrium in the cloud from a variety 
of initial abundances and intensities. Since we ignore the 
evolution of the background, there is some ambiguity to 
determining the minimum field necessary to drive II2 be- 
low detection limits in the cloud. Rather than address 
the many evolutionary possibilities for the background 
and whether it even arises given a set of initial circum- 
stances, we proceed from the assumption that it does 
appear. We then determine what intensities are required 
to suppress H2 formation below detection thresholds for 
a particular neutral H column density and metallicity. 
This stance is logical, given that the absence of II2 in 
the GRB sight lines renders the presence of LW radiation 
extremely likely and that its existence is of greater rele- 
vance than its origin, or history. However, note that in 
doing so we exclude some possibilities, such as large early 
H2 abundances capable of quenching LW backgrounds. 

We adopt a cloud core density Ucof 10 cm~^ and cloud 
densities rid of 5, 10, and 50 cm~'^. We examine metal- 
licities of 0.01, 0.1, and 1.0 Zq, for a total of nine models. 
Our coordinate grid was 200 uniform zones, with inner 
and outer boundaries of 0.1 and 1 kpc, respectively, and 
the same boundary conditions as in our ealier models. 
Our choice of mesh accelerates convergence to the H2 
abundances at some cost in accuracy to the H II region, 
but this is of little consequence to the molecular frac- 
tions in the cloud. In each case the simulation is run 
until H2 comes into equilibrium with the galactic back- 
ground. The H2 mass fraction is initialized to 2.0 x 
10^'*. Although not of direct importance to the molecu- 
lar abundances, the time-dependent spectrum of the 40 
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Mq cluster was used in these runs. 

The minimum background intensities required for H2 
suppression below the detection limit (defined to be 
fjjj =10~^-^) are plotted for all nine models in FigO It is 
clear that thresholds rise with both ambient density and 
metallicity as expected, varying from 0.144 to 144 times 
the Galactic UV mean. Initial H2 fractions greater than 
2.0 X 10~^ merely cause the cloud to reach equilibrium 
later and do not alter the final abundances. We have ver- 
ified that our results hold with initial fractions as large as 
0.2. The threshold intensities for 0.1 Zq < Z < Zq are 
consistent with the FUV backgrounds observed in star- 
forming regions in the LM C and SMC, whic h appear in 
the shaded region in Fig [5] ()Tumhnson et al.| [2002) . 

The scaling of threshold intensities with both cloud 
density and metallicity can be understood with simple 
equilibrium arguments. The timescale for H2 destruction 
in the cloud is the just the inverse of the photodissocia- 
tion rate coefficient 

idis = fcdis^\ (16) 

which is 1.84 x 10^^ s for the mean Milky Way UV back- 
ground. This is faster than any hydrodynamical process 
in the outer regions of the cloud, so final H2 fractions 
are determined by equilibrium between creation by dust 
(which completely dominates formation by the chan- 
nel at these densities) and destruction by LW photons. 
Taking maximum mass fractions f//^ = 10~^'^, one can 
approximate threshold intensities by equating the cre- 
ation and destruction processes: 

- 7.77x10^/h/ ^^^^ 

These equilibrium estimates, which we plot in FigO dif- 
fer from the code values by at most 30 %. The intensities 
scale linearly with metallicity and density, as confirmed 
by ZEUS-MP. 

A galactic LW background can be estimated for the 
afterglow spectra of TET07 using eq[T7l The density tih 
of the GRB host cloud can be fixed by assuming that 
the neutral H column density of the afterglow is local to 
the GRB and dividing it by some fiducial diameter for 
the cloud. We chose 10, 100, and 1000 pc for the size 
of the cloud and adopted Nh from column 3 of Table 1 
in TET07 for each spectrum. The detection limit f^a 
was computed from Njj and NH2,min from column 8 of 
Table 1. We plot the three intensity thresholds for II2 
destruction for 4 of the 5 GRBs in TET07 in Fig [5] as a 
function of metallicity (which are found in column 4 of 
Table 1). 

We find that the FUV backgrounds implied by eg [T71 
for the afterglows are consistent with those in our mod- 
els of comparable metallicity, particularly those associ- 
ated with more realistic diameters of 10 - 100 pc for the 
GRB host cloud. If we extrapolated our models to lower 
metallicity, they would also include the LW backgrounds 
associated with GRB 030323 for somewhat larger but 
still reasonable cloud sizes. Perhaps of most interest, 
afterglows with metallicities comparable to those of the 
SMC and LMC have similar uniform backgrounds, lend- 
ing weight to the notion that GRBs originate from ac- 
tive star-forming regions, even though their host galaxies 
cannot be observed. 



Our thresholds are valid for larger clouds at the same 
density due to the uniform nature of the background. 
The converse of these results is also true: if the host 
cloud has LW backgrounds below these levels, the GRB 
afterglow will exhibit H2 absorption regardless of the ini- 
tial molecular fractions. Since the borderline intensities 
are similar to those in regions of active star formation, 
at least in low metallicities, we anticipate it is only a 
matter of time before molecular hydrogen is detected in 
afterglows, if indeed GRB 060206 is not already an ex- 
ample of one (|Fvnbo et a l.| [2006) . Since we neglect both 
transport and self-shielding of LW flux, which reduce the 
dissociation of H2, our results somewhat underestimate 
the minimum backgrounds required to destroy molecular 
hydrogen. 

6. CONCLUSIONS 

The absence of H2 in the large neutral H columns of 
GRB afterglow spectra is easily explained by the H II re- 
gion of the progenitor or its cluster and by values of galac- 
tic LW backgrounds that are consistent with star-forming 
regions. In most cases the progenitor alone is sufficient 
to create the H II region. Since the H II region is illu- 
minated when the progenitor dies, and because stars in 
the cluster form at nearly the same time, the minimum 
radius at which neutral H is observed in the sight lines 
is determined by both the lifetime of the progenitor and 
sound speed in the ionized gas rather than by the number 
of stars in the group or its composition. As a result, these 
radii (as inferred from strong Mg I absorption) as a rule 
will be 50 - 100 pc. The sound speed (and therefore the 
expansion rate of the D-type front) is a relatively weak 
function of flux and spectral proflle, so we expect these 
radii to be a robust feature of afterglow spectra. The 
fairly narrow band in which they fall corresponds to the 
spread in MSL for likely GRB progenitors and is nearly 
identical for solitary stars and clusters, making it diffi- 
cult to discriminate whether the star was in a group or 
alone. Stars less massive than 25 Mq might fail to create 
an H II region, but they are not likely GRB candidates. 
Our findings are consistent with observed H II regions 
such as 30 Doradus in the LMC, who se radius is ~ 150 
pc an d contains around 50 OB stars (|Massev fc Hunted 
Il998h . 

Total flux plays a role in the sizes of the H II regions 
only if it is so intense that the front is still R-type at the 
time of the burst. Unrealistic numbers of massive stars 
would be necessary for this to happen, especially if dust is 
present. Winds, although not included in these models, 
will reinforce our range of minimum radii by reducing 
their dependence on the density of the cloud core. X- 
rays and extreme UV from the wind-blown shell may 
alter the structure of H II regions but not their radii or 
II2 abundances in the cloud. 

H2 abundances in the ambient ISM are nearly indepen- 
dent of the properties of the GRB progenitor or host clus- 
ter. They are set by the mean intensity of the UV back- 
ground, dust abundance and gas density. Since the min- 
imum intensities required to extinguish II2 are common 
to regions of vigorous star formation, it is not surprising 
that molecular hydrogen has not yet been conclusively 
identified in GRB afterglows. On the other hand, these 
backgrounds are large enough that some clouds might 
not achieve them, so it is likely that II2 will eventually 
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be found along a GRB sight line. 

All of our models hold the GRB progenitor at the cen- 
ter of the molecular cloud core, but this in reality is un- 
likely due to the peculiar motions of the stars in the clus- 
ter. However, stellar motion within the cluster would be 
at most 10 pc and would still be bracketed by our models 
because the progenitor remains deeply imbedded within 
the cloud. If the star wandered toward the edge of the 
core, the size of its H II region would depend only on 
rici, not ric. As our models show, the radius would not 
increase by more than 10 % in these circumstances. The 
I-front would blow out the wall of the core, spilling ion- 
ized gas o ut into the c loud in a supersonic 'champagne' 
flow {e.g. lYorke et al.| [T984l. 

One issue that remains unclear is the effect of the op- 
tical/UV afterglow itself on molecular hydrogen in its 
path. Alt hough momentary, it i s extremely intense (10^° 
photons). iDraine fc Had (f200l) convincingly show that 
the optical flash in neutral molecular cores will destroy 
H2 only out to 5 - 10 pc and have little effect on molecu- 
lar fractions in the rest of the host cloud. However, this 
treatment does not account for the H II region, which 
hollows out a diffuse cavity of gas ~ 100 pc in radius. 
What remains to be determined is if the optical pulse 
can destroy any H2 beyond the H II region that was not 
dissociated prior to the burst. If so, this would lower the 
minimum LW background necessary to extinguish molec- 
ular hydrogen in the rest of the cloud; if not, our results 
would unambiguously prove that the H2 was absent be- 
fore the explosion. This question cannot be settled by 
the current version of our code because the static ap- 
proximation to radiative transfer we employ breaks down 
in the limit of these short, intense flashes. Upgrades to 
our algorithm that retain the time-dependent term of the 
transfer equation will enable us to resolve these issues in 
a future study. 

Multidimensional simulations now in progress will re- 
veal whatever dynamical instabilities occur in the ion- 
ization front and wind bubble. Efficient radiative cool- 
ing by both metals and H2 in the shocked neutral gas 



accumulated by the I-front wi ll collapse it into a thin 
shell prone to fragmentation (fOarcia-Segura fc Francol 
ll996HWhalen fc Normanll2007al) . Ravleigh-Tavlor insta- 
bilities may appear in the wind bubble even in the ab- 
sence of cooling if it accele rates in density gradients [e.g. 
iMac Low fc Normanlll993l ). Both phenomena can clump 
gas on multiple spatial scales in the flow prior to the 
burst and might be imprinted on afterglow spectra. Fu- 
ture observations may uncover these features and yield 
insights into the environment of the burst. 

Finally, it will be valuable to compare the incidence 
of H2 detections in GRB host galaxies with the inci- 
dence measured in QSO damped Lyman-alpha (DLA) 
sight lines. The latter show an inci dence of ~ 10% for 
molecular fractions /(H2) > 10~'*'^ (|Noterdae mc et aTl 
[20 08 ). Current expectation is that the QSO-DLA sight- 
lines penetrate galaxies that hav e lower total SFRs at a 
larger average imp act parameter (jProchaska et al. II2007I : 
iFynbo et al.ll2008h . The results presented here suggest 
that one may test these hypotheses, in part, by mea- 
suring the relative incidence of H2. For example, one 
may expect a lower average UV intensity along QSO- 
DLA sightlines than GRB-DLA sightlines and therefore 
a higher incidence of H2 detections. 
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